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1.  INTRODUCTION 

When  a  center-driven  dipole  antenna  is  placed  parallel  to  a  conducting 
or  dielectric  half-space  at  a  height  d  above  it,  distributions  of  current  and 
charge  are  induced  in  it  that  interact  with  the  currents  and  charges  in  the 
antenna  and  combine  with  these  to  generate  the  electromagnetic  field.  The 
calculation  of  the  interaction  is  difficult  except  when  the  half-space  is 
perfectly  conducting.  In  this  case  the  effect  of  the  distributions  of  cur¬ 
rent  and  charge  on  the  conducting  surface  is  precisely  that  of  the  equal  and 
opposite  currents  and  charges  in  a  fictitious  “image"  antenna  located  at  a 
distance  d  below  the  surface.  It  replaces  the  entire  half-space  and  reduces 
the  determination  of  currents,  charges  and  electromagnetic  fields  to  that  of 
two  identical  parallel  antennas  separated  by  a  distance  2d  and  driven  by 
equal  and  opposite  generators.  When  2kgd  <<  1,  the  antenna  and  its  image 
constitute  a  two  wire  transmission  line  with  currents  and  driving-point  ad¬ 
mittance  that  are  determined  primarily  by  the  impedance  and  capacitance  per 
unit  length  of  the  line  and  affected  only  negligibly  by  radiation.  This  is 
true  of  the  several  circuits  shown  in  Fig.  1:  the  center-driven  line  with 
open  ends  shown  in  Fig.  la;  the  line  terminated  at  each  end  In  an  arbitrary 
shunt  impedance  ZL/2,  where  7^  may  be  the  characteristic  impedance  of  the 
line,  as  in  Fig.  lb;  or  the  line  with  a  series  impedance  7^/2  at  a  quarter 
wavelength  from  the  open  end  as  in  Fig.  lc.  In  each  case  the  image  conduc¬ 
tors  and  impedances  are  shown  in  broker,  lines.  Note  that  the  circuits  of 
Figs,  lb  and  lc  are  equivalent  insofar  as  the  current  and  driving-point  ad¬ 
mittance  of  the  main  line  between  the  impedances  Z^/2  are  concerned. 

When  the  half-space  below  the  center-driven  dipole  Is  not  perfectly  con¬ 
ducting  but  consists  of  a  material  like  the  earth,  the  sea  or  a  fresh  water 
lake,  conduction  or  polarization  currents  or  both  of  these  may  be  induced. 
However,  they  are  not  confined  to  a  thin  surface  layer  but  penetrate  more  or 
less  deeply  into  the  medium.  Their  effect  on  the  current  in  the  antenna  and 
their  contribution  to  the  elec tremagne tic  field  are  not  easily  determined. 

In  the  treatment  of  antennas  parallel  to  the  earth  it  has  been  common  prac¬ 
tice  to  assume  the  distribution  of  current  to  be  known  and  to  be  substantial¬ 
ly  that  for  the  same  antennas  when  isolated  (1J.  This  implicitly  also  takes 
for  granted  that  the  wave  number  for  the  current  on  the  highly  conducting  an¬ 
tenna  is  the  same  as  that  for  the  isolated  antenna,  which  is  the  free-space 
wave  number  kg  *  w/c  where  c  is  the  velocity  of  lip  it.  Perhaps  the  most  im- 
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porfcant  antenna  of  this  type  is  the  Beverage  antenna  which  consists  of  a  wire 
parallel  to  and  dote  to  the  earth  terminated  in  a  manner  to  provide  a  trav¬ 
eling  wave  of  current.  This  antenna  was  first  described  by  Beverage  [ 2 J .  It 
has  been  discussed  by  numerous  authors  primarily  as  a  directional  receiving 
antenna  [3]~[6J.  In  general,  a  oimple  traveling  wave  with  the  free-spaee 
wave  number  has  been  assumed  and  the  dependence  of  the  wave  number  and  of  the 
distribution  of  current  on  the  properties  of  the  adjacent  half-space  has  not 
been  determined.  Wait  [7]  derived  a  modal  equation  for  the  wave  number  for 
an  assumed  traveling  wave  along  an  infinitely  long  wire  parallel  to  a  dissi¬ 
pative  half-space,  but  concluded  that  its  solution  in  the  general  case  is  a 
"formidable  task"  and  that  "in  spite  of  the  fundamental  nature  of  the  problem, 
there  has  not  been  a  great  deal  of  progress  in  obtaining  useful  results." 

Even  in  the  limiting  transmission-line  form.  Wait  expresses  his  final  formula 
for  the  wave  number  in  terms  r»f  an  unevaluated  integral. 

2.  REVIEW  OF  SOLUTION  FOR  CURRENT  IN  ECCENTRICALLY  INSULATED  ANTENNAS 

In  a  recent  paper  [8]  a  new  solution  was  derived  for  the  current  in  and 
the  admittance  of  a  center-driven  insulated  dipole  when  immersed  in  an  ambi¬ 
ent  medium  with  a  characteristic  complex  wave  number  much  greater  in  magni¬ 
tude  than  that  of  the  insulating  material.  This  theory  was  verified  by  very 
extensive  measurements  of  distributions  of  current  and  charge  along  air- 
insulated  antennas  immersed  in  a  large  lake  [9].  The  theory  of  the  insulated 
antenna  has  now  been  generalized  to  include  eccentric  insulators  [10].  That 
Is,  the  antenna  is  not  located  along  the  axis  of  the  circular  cylinder  of  in¬ 
sulation  but  is  displaced  by  an  arbitrary  distance  D.  The  cross  section  of 
an  eccentrically  insulated  antenna  is  shown  in  Fig.  2.  The  radius  of  the  in¬ 
sulation  is  b,  that  of  the  conductor  is  a.  The  axis  of  the  conductor  is  at  a 
distance  D  from  the  axis  of  the  insulator  or  a  distance  d  *  b  -  D  from  its 
surface.  The  wave  number  characteristic  of  the  insulating  material  in  the 
range  a  <  r  <  b  is  +  i'^  *  +  io^/w) j  7 *"  *  m(p i  the  wave 

number  characteristic  of  the  ambient  medium  in  the  range  b  <  r  <  «  is  k,  * 

1/2  -  1/2  H 
8^  +  ia^  ■  aj[p(e^  +  io^/aOj  »  where  6^  is  the  real  phase  con¬ 

stant,  a,  the  real  attenuation  constant.  It  is  assumed  in  the  analysis  that 
the  following  two  conditions  are  satisfied: 

(k2d|  «  1 


i  k4 i  >>  |k2 1 
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FIG.  2  CROSS  SECTION  OF  ECCENTRICALLY  INSULATED 
ANTENNA 
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Subject  to  these  conditions  the  distribution  of  current  along  the  eccentri¬ 
cally  insulated  dipole  of  half  length  h  is 

,e 


I(z) 
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22  cos  k.  h 
c  L 


(2) 


where  a  tine  dependence  exp(-iwt)  has  been  assumed  and 
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The  complex  wave  impedance  of  region  2  is  ^  *  (v/ty)  *  1°  the  last  ex¬ 

pression  , 
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This  locates  an  infinitely  thin  wire  as  shown  in  Fig.  2. 
The  driving-point  admittance  is 


Y  *  G  -  iB 


1(0) 


27 


tan 


h" 


(8) 


Thi3  is  the  admittance  of  two  open-ended  sections  of  transmission  line  with 
characteristic  impedance  Z  and  complex  wave  number  kf  in  series  with  two 

C  L. 

generators  each  with  driving  voltage  V^/2  at  z  ■  0  as  Indicated  in  Fig.  3. 
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3.  THE  COMPLEX  WAVE  NUMBER  FOR  THE  ANTENNA  NEAR  A  PLANE  EARTH 

The  wave  number  for  the  antenna  eccentrically  placed  in  a  circular  cyl¬ 
inder  of  radius  b  is  given  by  (4) .  In  order  to  obtain  its  value  in  the  limit 
b  -*•  D  «  with  d  «  b  -  D  fixed  and  finite,  it  is  necessary  to  determine 

and  The  former  is  readily  done.  Thus, 


cosh 
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2  2 

where  the  logarithmic  approximation  is  valid  when  a  <<  d  .  The  evaluation 
of  Aq  is  more  difficult.  As  a  first  step,  note  that  with  (7) 
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so  that  (6)  becomes 
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or,  in  a  form  more  convenient  for  later  use. 


A  «  lim 
u  b-«“ 


lim 

b***0 

D-w» 


H£1>(k*b) ,  j*-*y2  ?  ( 

{b  -  d 

\2m  H  , (k.b) 
j  m-1  4 

Hll,(k4b)  \  b  /  m*2 

V  b 

/  H(1>(k,b) 

m  4 

00 

i  +  2  I  1 

m-2 

'b  -  d 

\2m  H(1}(k.b) 

\  m-l  4 

V  b 

1  H*15 (kb) 

m  4 

(12) 


It  is  not  permissible  to  take  the  limit  b  m  before  the  evaluation  of  the 
sum,  since  this  involves  m  ■*  This  means  that  the  familiar  asymptotic 
forms  for  the  Hankel  functions  of  large  argument  cannot  be  used  for  m  >  Ik^bj 
since  they  are  derived  on  the  assumption  that  the  order  is  small  compared  to 
the  argument.  They  have  been  used  for  orders  0  and  1.  What  is  required  is 
an  expression  for  the  Hankel  functions  that  is  valid  for  both  large  orders 
and  large  arguments.  Such  a  form  may  be  obtained  with  the  recurrence  rela¬ 
tion 
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zH(i:>,(z)  +  mH(1)(z)  «  zH(1J(z) 
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which  yields 
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where  the  prime  denotes  the  derivative  with  respect  to  the  argument.  The 
following  results  are  obtained  with  formulas  given  in  Bateman  [11,  p.  86] : 
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(15) 
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When  the  differentiation  is  carried  out. 
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In  each  of  these  expressions  the  first  term  is  neglected  in  the  final  form 
since  it  is  small  when  the  inequalities  on  the  left  are  considered.  The  sub¬ 
stitution  of  these  results  in  (14)  yields: 
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With  (19)  and  (20),  (12)  becomes: 
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The  sums  in  this  expression  can  be  converted  into  integrals  with  the  substi¬ 
tutions  X  ■  m/k^b,  dx  *  1/k^b.  Also, 
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With  these  substitutions,  (21)  is  approximated  by: 
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where  A  =  2k4d.  This  expression  can  be  rearranged  as  follows: 
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The  several  integrals  in  (26)  are  evaluated  in  the  Appendix.  The  final  ex¬ 
pression  for  Ag  is  obtained  with  (A-17).  It  is 
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where  A  *  2k, d,  with  d  the  distance  from  the  axis  of  the  «r;.>.enna  in  air  to 
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the  plane  boundary  between  the  air  and  the  material  half-space.  Note  that  in 
the  limit  b  ■*  »,  the  first  term  in  (26),  i,  vanishes. 

The  substitution  of  (5)  and  (27)  into  (A)  gives 
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Thus,  the  complex  wave  number  k  for  the  antenna  in  air  over  a  plane  medium 

L 

with  characteristic  wave  number  k^  depends  only  on  the  distance  d  from  the 
axis  of  the  antenna  to  the  plane  boundary,  the  radius  a  of  the  antenna  and 
the  wave  numbers  k^  and  k£  “  k^  of  the  air.  Note  that  if  j  2k^d |  <  1,  only 
two  terras  in  the  series  in  (28)  are  required.  For  example,  when  k^  is  real 
and  k^d  *=  0.5,  1^(1)  -  0.6019,  I  (1)  -  0.5652  and 
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With  the  complex  wave  number  k^  determined,  the  characteristic  impedance 
is  obtained  from  (3)  with  (9)  and  (28).  Thus, 
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The  distribution  of  current  along  the  antenna  is  given  by  (2),  the  driving- 
point  admittance  by  (8),  with  k^  and  Z£  defined  in  (28)  and  (30). 

The  effective  line  constants  per  unit  length  are  readily  evaluated  since 
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where  and  are  the  real  effective  conductivity  and  permittivity  of  the 
insulating  cylinder.  Similarly, 
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usv, 


with 


e  i  e  i 

z_  *  z  +  z.  *  -iul  +  z. 
LA  4 


,e  p  _  2d 
i  *  in  — 
2v  a 


i  -iup 


{_  (2k4d)2  2V 


K1(2k4d)  ini 


l(2k4d)  / 

Ak4d  "  V 


(34) 


(35) 


2k, d  (2k, d) 3  (2k4d)5 

3  +  45  *  1575  +  ** 


06) 


Thus,  the  antenna  over  a  half-space  characterized  by  |k4|  >>  | k^ !  behaves 
like  a  transmission  line  with  the  line  constants  given  in  (32),  (35)  and  (36), 
If  the  ohmic  losses  in  the  antenna  are  to  be  included,  an  additional  imped¬ 
ance  per  unit  length,  viz,, 


i  V-W0! 


2na 


(37) 


must  be  added  to  z^  in  (34)  and  included  in  k^  *  -z^y^  an<^  2c  *  ^ZL^L’  *n 
(37),  P^  and 
the  antenna. 


(37),  P^  and  are  the  permeability  and  conductivity  of  the  metal  forming 


4.  THE  AXIAL  ELECTRIC  FIELD  ON  THE  SURFACE  OF  THE  HALF-SPACE 

The  transform  of  the  axial  component  of  the  electric  field  at  a  radius 
r  from  the  axis  of  the  insulating  ragion  2  that  contains  the  eccentrically 
located  antenna  is  [10] : 


E  (r,6)  - 
Z  *  2*U£ 


(k2  -  C2)n(r,6)  +  r-r-  A(r,6) 

2  L  4° 


(38) 
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vhere  0  It  the  angle  aeasured  from  the  radius  through  the  antenna  where  this 
is  closest  to  the  boundary,  I  is  the  transform  of  the  current,  and 


..  +  ,  r  (  "o  \"  ;*rx_b) 

“  [u  +2  2.  I  2  J  m  \ 

Hj1J(k4b)  m«l\  b^  /  H^J(k4b 

,  f  (rxn)2  -  2b2rxn  cos  8  +  b4  1 
fi(r,6)  -i*n  ^ - j- 

{_  b  (r  -  2xQr  cos  6  +  xQ)  J 


cos  m6 


At  the  boundary  r  »  b  and  with  (10),  Xq  -*■  b  -  d  as  b  -*■  Hence, 
fl(b,0)  -  0 


Mb,  6) 


H0U<V>  +2  J  / 

'b-d^ 

m  H(1) 

^^(k.b)  jl' 

k  b  J 

1  H(1> 

m+1 

<V»> 

(k,b) 


cos  mO 


It  follows  that 


Ez(b,0) 


±  iwpl  A(b,0) 

2  it  k.b 

A 


which  is  independent  of  the  transform  variable  c.  Hence,  the  inverse  Fourier 
transform  with  respect  to  z  gives  simply: 


Ez(b,0) 


iupl  A(b,0) 


for  the  z-component  of  the  electric  field  in  terms  of  the  current  I. 

The  evaluation  of  A(b,0)/k4b  can  be  cart-led  out  by  the  procedure  used 
for  A^/k^b.  With  the  notation  L  *  b0  (where  L  is  the  distance  along  the 
boundary  from  0  *  0  in  directions  at  right  angles  to  the  antenna),  B  *  k^L 
and  A  «  k^d  (not  2k^d  as  in  the  evaluation  of  A^/k^b),  the  equation  corre¬ 
sponding  to  (25)  is 


A(b,0)  *  i  +  2k -b j  J  (x  +  iVl  -  x“)e*^x  cos  Bx  dx  +  /  (x  -  Vx^  -  l)e 

L  o  1 


*  cos  Bx  dx 
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(Note  that  this  reduces  to  (25)  when  B  -  0.)  The  evaluation  of  the  square 
bracket  in  (45)  Is  cert led  out  in  the  Appendix,  and  given  in  its  general  form 
in  (a-13) .  In  the  limit  b  «,  the  first  tern  in  (45)  vanishes.  The  final 
result  is 


lim 

b-K» 


Mb,  8) 

V 


D-«>> 


A(d,L) 

k4b 


A2  -  B2  K1(A  "  1B)  .  _* 

; '2  ~~272 - A  -  iB  +  ““  2(B  +  iA) 

(A  +  B  ; 

X  tEx(B  +  iA)  +  Y1(B  +  iA)  ~ 


(46) 


where  Kj^  is  the  modified  Bessel  function  of  the  second  kind  and  order  one, 
is  the  Weber  function  of  order  one  and  Y1  is  the  Bessel  function  of  the 
second  kind  and  order  one.  These  are  defined  in  the  Appendix.  Accordingly, 
with  A  »  k^d,  B  •  k^L, 


E  (d,L) 
z 


itoyl 


a2-b2 

\  (A2  +  B2)2 


KX(A  -  iB) 


-  -  A -IF-  +  iItt  2(bTW  [ei(b  +  1A) 


+  VB+  ^ 


(47) 


At  L  *  0,  (A-17)  gives: 


Er(d,0) 


K.(A)  in I. (A) 

— - j. - i - 

A  2A 


(48) 


This  is  the  axial  field  at  the  boundary  surface  along  the  normal  projection 
of  the  antenna  as  shown  in  Fig.  4. 

Since  these  expressions  are  not  particularly  transparent,  it  is  advan¬ 
tageous  to  examine  them  it.  certain  limiting  cases.  Consider  first  their 
forms  for  large  arguments,  i.e.,  when  |A  +  B  j  >>  1.  The  applicable  formula 
is  (A-19).  Thus, 


E2(d,L) 


iwpl 

A2  -  B2  .1  .  ^ 

f  *  -(A-  iB) 

[l+  3 

IT 

((a2  +  b2)2  A-1B* 

2 (A  -  iB)  6 

[/  +  8 (A  -  iB) 

15 


128 (A  -  iB) 


+  .  • , 


-  lira 


B  +  IA 


1  + 


(B  +  iA)  (B  +  iA) 


(coot.) 


ka 


k4 


2a 

-*i  CONDUCTOR  (1) 

AIR  (2) 

L 

- •»£ - 

Ez (down) 

AMBIENT  MEDIUM  (4) 


FIG.  4  HORIZONTAL  CONDUCTOR  (1)  IN  AIR  (2)  OVER 
AMBIENT  MEDIUM  (4);  Ez  AT  BOUNDARY 
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45 


(49) 


(B  4-  iA) 

/ 

the  field  at  L  •>  0  for  |a|  »  1  is  obtained  with  (A-20).  It  is 


E2(d,0) 


iupl 


\±-l  JZe“A 

1  ,2  A  I2A 


u 


1  +  5X“ 


15 


12  8A 


***J  +  a[x  “ 


(50) 

A  A  9  9  9 

The  condition  {Az  +  B  |  *  |k^(L  +  d  )|  >>  1  can  be  satisfied  in  several  ways. 
When  Jk^|  is  made  very  large  -  or  infinite  for  a  perfect  conductor  -  the 
field  along  the  surface  of  the  half-space  is  very  small.  For  ■* 


Ez(d,L)  -  0 


(51) 


as  it  should.  When  ( A {  *  |k^d|  is  very  large  because  |k^|  is  large,  as  for  a 
good  conductor,  Ez(d,L)  is  very  small.  For  a  given,  sufficiently  large  k^d, 

the  electric  field  along  the  boundary  decreases  with  increasing  B  *  k.L  ac- 

2  2  4 
cording  to  (49).  When  B  >>  A  ,  the  field  is  obtained  with  (A-21).  It  is 


*  ittlil  (  1  i  JlH  iB  I”.  31 

:(d*L)  “  12&  e  [1  +  SB- 


15 


128B 


1) 


(52) 


Thusl  for  large  values  of  B  *  k,L,  the  field  decays  radially  outward  along 

1/7 

the  boundary  surface  as  l/(k^L)  '  . 

2  2  2  2  2 

The  condition  for  small  arguments  |k^(L  +  d  )|  «  |A  +  B  |  «  1  is 

quite  restrictive  but  not  physically  unavailable.  Note  that  the  conditions 

|k^|  >>  (k^l  and  (k^d j  <<  1  have  been  imposed  from  the  outset.  To  require, 

in  addition,  that  |k2d2|  «  1  implies  that  jk^dj  <<<  1.  Specifically,  if  k^ 

applies  to  lake  water  with  ■  81,  it  is  given  by  k^  *  9w/c  *  3w  *  10“®  m~^. 

If  the  antenna  is  located  at  a  distance  of  d  “  1  m  from  the  surface  of  the 

—8 

water,  k,d  *  3w  x  10  .  For  this  to  be  small  compared  with  one,  viz., 

|k,d|  <  0.1,  a)  must  not  exceed  about  3  x  10°.  Note  that  k-d  •  u>d/c  <  0.01, 

**  “  1/2  ^  ~ 
Similarly,  for  sea  water  with  o  »  4  mho/m,  k^  -  (iupo)  and  |k^|  » 

4(uir  *  10~7)1/2.  With  d  *  1  m,  |k^d|  «  4(u>ir  *  10~7)*72.  This  will  satisfy 
the  condition  jk^d|  <  0.1  if  u  does  not  exceed  about  2  x  10^.  At  this  fre¬ 
quency  kjd  -  6.67  x  10”6,  k^d  *  0.1eilt7A. 

2  2 

Subject  to  the  condition  |A  +  B  j  «  1,  an  approximate  formula  for 
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E  (d,L)  is  obtained  with  (A-27).  If  is 

Z 


E2(d,L)  * 


iupl 

2ir 


(53) 


It  is  interesting  to  compare  this  expression  with  the  incident  field  gener¬ 
ated  by  the  current  in  the  wire  if  this  is  treated  as  a  line  source.  It  is 

E^nC(d,L)  -  iu pi  ~  +  d2)  (54) 

2  2  2 

Sufficiently  near  the  antenna  where  k2(L  +  d  )  <<  1,  the  snail-argument  form 
of  the  Hankel  function  is  a  satisfactory  approximation.  This  gives: 

-  Y  +  1  j  (55) 


The  total  field  E  (d,L)  at  the  boundary  as  given  in  (53)  can  be  ex- 
z 

pressed  as  the  sum  of  the  Incident  field  plus  the  field  scattered  from  the 
matter-filled  half-space.  With  (53)  and  (55)  it  is 


Ez(d,L)  -  EznC(d,L)  +  E*CS£ 


I 


I L2  f  d2 


Y  +  1 1  ” £n  iq  + 1 


This  formula  Indicates  that  at  least  within  the  range  of  the  conditions 
2  2 

|k^j  >>  Jkjl  and  |k^d|  <<1,  the  contribution  of  the  currents  and  charges 
in  the  matter-filled  half-space  to  Ez(d,L)  on  its  surface  is  the  scattered 
field 


(57) 


This  field  cannot  be  interpreted  as  being  due  to  an  imaga  current  in  an 

image  antenna  at  a  distance  2d  from  the  antenna  since  it  is  independent  of  d. 

£ 

For  lake  water  with  er  *  80,  d  ■  1  n,  w  “  3  *  10  ,  k^d  ■  0.01,  k^d  * 
0.09,  the  ratio  of  the  scattered  to  the  Incident  field  is: 


(58) 


Escat(dtO)  fcn(k4/k2)  -  0.5  12.82 

Einc(d>0)  “  "  ln(2/k2d)  -  Y  iir/2  “ 


0.35e 


il&2* 


The  ratio  of  the  total  field  to  the  incident  field  is: 


S&OL.  .  0.68eiO-16  «  0.68e19*-5 

EinC(d,0)  EinC(d,0) 


(59) 


V 


For  sea  water  with  o  »  4  mho/m,  d  *  1  m,  w 
*  O.le**7^,  so  that 


2  x-  io k2d 


6.67  x  10 


-6 


Escat 

Einc 


(d,0) 

(d,0) 


0.75e 


13.09 


0.75eil77e 


(60) 


and 


■  .  0.25ei0*16 

Einc(d,0) 


C.25e 
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(61) 


It  is  seen  that  the  effect  of  the  dielectric  (lake  water)  like  that  of  the 
conductor  (sea  water)  is  to  combine  with  the  incident  field  to  produce  a 
total  field  at  the  boundary  that  is  substantially  smaller  chan  the  incident 
field.  With  lake  water,  the  incident  and  scattered  fields  differ  in  phase  by 
162*;  with  sea  water  the  phase  difference  is  about  177*  for  the  particular 
distance  d  and  frequency. 


5.  THE  GENERALIZED  BEVERAGE  ANTENNA;  CONCLUSION 

The  well-known  Beverage  antenna  [1J~[4]  consists  of  a  wire  up  to  several 
wavelengths  long  that  is  stretched  parallel  to  and  an  electrically  short  dis¬ 
tance  from  the  surface  of  the  earth  to  which  it  is  connected  at  both  ends 
through  suitable  resistive  terminations.  Such  an  arrangement  to  achieve  a 
traveling  wave  is  possible  only  when  the  earth  is  sufficiently  conducting  to 
provide  an  adequate  ground  connection.  A  more  generally  useful  method  of  ob¬ 
taining  a  traveling  wave  along  a  sufficiently  long  wire  stretched  parallel  to 
a  material  half-space  of  any  material  with  properties  that  range  from  those 
of  a  good  dielectric  (lake  water)  to  those  of  a  good  conductor  (sea  water)  is 
to  connect  the  resistive  terminations  in  series  with  the  wire  at  distances 
near  a  quarter  wavelength  from  the  ends  as  in  Fig.  lc.  An  analysis  of  this 
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generalized  form  of  the  Beverage  antenna  in  the  framework  of  the  present 
theory  is  reserved  for  another  report. 

In  this  report  a  new  approach  has  been  given  for  the  traveling  wave  an¬ 
tenna  over  an  imperfectly  conducting  or  dielectric  half-space.  The  complex 
wave  number,  the  distribution  of  current,  and  the  electric  field  tangent  to 
the  boundary  have  been  determined.  The  problem  of  calculating  the  radiation 
field  of  such  an  antenna  with  the  known  distribution  of  current  can  be  accom¬ 
plished  using  well  known  methods. 
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APPENDIX:  EVALUATION  OF  INTEGRALS 

The  following  integrals  occur  in  (26)  with  B  ■  0  and  A  *  2k^d  and  in 
(45)  with  B  "  k^L  and  A  *  k^d : 


I  *  /  (x  +  i Ml  -  x2)e~Ax  cos  Bx  dx  +  / 
0  1 


CD  - 

+  /  (x  -  ?x2  -  l)e~^x  cos  Bx  dx  (A-l) 


This  expression  can  be  rearranged  as  follows: 

I  -  7  xe-Ax  cos  Bx  dx  +  j\ /  iVl  -  x2  e“(A~  1B)x  dx  +  /  (  Vx2  -  1  e“(A"  iB)x  dx 

0  zLo  1 


-7  Vx2-le-<A-1B>x  dx+i[/ 

i  Ho 

-/  Vx2  -  1  e-(A+iB)x  dx] 


.  2  -(A+  iB)x  , 

i  TX  -  x  e  dx 


(A-2) 


Here  the  first  integral  is: 

®  2  2 

1  «  /  xe^*  cos  Bx  dx  -  -£-  ~.  A  (A-X 

0  (A^  +  bV 

With  the  substitution  x  *  cosh  6,  the  fourth  integral  is  evaluated  as  follo\  j, 

_____  JB\V  O  _/*  _  i  ON  a  ^1  (A  —  iB) 


I,  -  /  Vx2  -  1  e-(A'iB>x  dx  -  /  sinh2  6  e-<A'*B)cosh  6  dfi  g 
1  1  0  A  -  IB 

(A-4) 

where  :.^(z)  is  the  modified  Bessel  function  of  the  first  order  and  second 
kind.  Hie  remaining  integrals  are  equivalent  to 

I,  -  if  /  iVT77  d„  +  /  ViXTT  e-(A-  “>■  dx]  -  [  iVTX? 

Lo  i  J  [o 


-(A-iB)x  .  .7  VT77„-(A-iB)x 


dx  +  j  fx  -  1  e" 
1 


(A-5) 


where  the  asterisk  denotes  the  complex  conjugate  of  the  expression  in  the 
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brackets,  Since  the  second  square  bracket  in  (A-5)  is  the  complex  conjugate 
of  the  first,  their  difference  is  twice  the  imaginary  part  of  the  first. 

Thst  is. 


I*  -  lira  /  i  Vl  -  x2  e“(A"lB)x  d x  +  /  V>?  -  1  e" 
J  0  1 


(A  -  iB)x 


(A-6) 


This  is  an  integral  along  the  real  axis  from  zero  to  infinity  with  an  upward 
deflection  around  the  point  x  *  1  and  the  attached  downward  branch  cut.  The 
path  of  integration  can  be  rotated  counterclockwise  90°  to  the  y-axis  with 
the  substitution  x  *  iy.  This  leaves  simply 

I  «  ilm  /  i  Vl  +  y2  e“(A“  ldy  »  _iIm  J  cosh2  e  e“(B  +  6  de  (A_7) 

J  0  0 

The  last  integral  in  (A-7)  is  readily  identified  as  the  sum  of  the  following 
integrals  with  suitable  coefficient: 


Vz)  +  Vz>  -  - 1  /  * 


2  r  -z  sinh  9 


(A-8) 


E2(z)  +  Y2(z) 


/  cosh  20  e 


-z  sinh  0 


(A-9) 


E  (z)  is  the  Weber  function  of  order  n  [31,  p.  35],  [12],  Y  (z)  is  the  Bessel 
n  n 

function  of  the  second  kind  and  order  n.  With  the  well-known  recurrence  re¬ 
lations 


Yq(z)  +  Y2(z)  -  f  Y 1(z) 


(A-10) 


E0U>  +  e2(z)  "  f  [El<z>  'I] 


(A-ll) 


(A-7)  becomes 


h  ■  1Im  hbttst[ei(b  -  “>  +  V”  +  “>  -  !] 


(A-12) 


When  (A-3),  (A-4)  and  (A*1 2)  are  combined  to  give  (A-2) ,  this  becomes 
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1  *  X1  -  X2  +  *3 


I  - 


A2  -  B2 
(A2  +  B2)2 


KX(A  -  iB) 

A  -  IB  +  1Im  T(B* 


[  E1(B  +  1A)  +  Y1(B  +  1A) 


-I] 


(A-13) 


Special  Case  with  B  ■  0: 

An  important  special  case  Is  when  B  *  0,  In  this  case 


K  (A)  itrl  (A) 

k2  A  *  2A 
A 


-fs[El(iA>  -I] 


(A-14) 


where  I^(A)  is  the  modified  Bessel  function  of  the  first  kind  and  first  order. 
This  follows  with 


Y.CiA)  /j\  1  r  2  ~\ 

Im  ■  -  Im  »  [-Hji;(iA)  +  J^iA)]  -  Im  K^A)  +  iI1(A)j  -  (A-15) 

and  [11,  p.  35] 

*!<*>  -  I 

Note  that  K^(A)  is  real.  Since  the  modified  Bessel  functions  of  real  argu¬ 
ments  are  available  in  Tables,  a  convenient  form  of  (A-14)  is: 


2  4 

1  3  45 


z 


8 


1575  99225 


(A-16) 


,  K.  (A)  iitl  (A) 

T  m  _  £  .  ...  +  _ = - i 

.2  A  2A  1 

A 


A  A3 
3  *45 


1575  99225 


+  .  •  < 


(A-17) 


It  is  this  form  that  is  used  to  obtain  (27)  from  (26) . 

2  2 

Special  Case  for  Large  Arguments;  A  +  B  >>1: 

When  the  argument  is  large,  the  following  series  is  useful  [12]: 


E1(z)  +  Y1(z)  % 


(A-18) 


With  (A-18)  and  the  usual  expressions  for  Bessel  functions  with  large  argu¬ 
ment,  it  follows  that 
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.2  „2 


I  -  -  A  -  "  B-  _  1  J - 7 

(A2  +  B2)2  A  “  iB  “ 


W 


-(A-iB) 


i  4. _ 3_  15 

8(A  -  IB)  - 


+  . . 


128 (A  -  IB)4 


-  ilm 


— — —  [  i  + _ i _ 3  _ 45 _  1 

L  (B  +  1A)2  (B  +  iA)4  (B  +  1A)6  *  J 

When  |a|  »  1,  B  «  0, 


(A-19) 


I  - 


A  L  “*  128A“ 

On  the  other  hand.  If  jB|  »  1  and  |b|  »  |a| , 

1  -I  Ijr- 

I  * 


JL  1. 

A2  ’a*"- 


J  (A-20) 


-  X  -  *  J*L  e1B  h  +  3i  15  1 

b2  bV2B*  L 1  +  ®B  +  i28B^  +  *  *  *  J 


(A-21) 


Special  Case  for  Small  Arguments;  jB2  +  A2j  <<  1; 
With  snail  arguments  [13] 


V*)  'f*nf  +  i  +  J(2Y- 


1) 


(A-22) 


(A-23) 


*V2)  ^  2(y  +  in  |  | 

so  that  with  (A-16) 

2T['l<!)  +E1(j>  + 

The  use  of  <A-21)-(A-23)  in  (A-13)  gives: 

-  J  [*”  +  V  -  y  +  M»[t(r  +  to  -  i  ) 


(A-24) 


I  ^ 


(A2  +  B2)2  (A 


- _ 1  _  B  +  iA  1 

(B  +  iA)2  3  J 


(A-25) 


This  reduces  to: 
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Here  the  leading  real  and  imaginary  terms  are 
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